Hybrid gap structure in the heavy-fermion superconductor Celrln 5 
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The thermal conductivity k of the heavy-fermion superconductor Celrlns was measured as a 
function of temperature down to T c /8, for current directions perpendicular (J || a) and parallel 
(J || c) to the tetragonal c axis. For J || a, a sizable residual linear term Ko/T is observed, as 
previously, which confirms the presence of line nodes in the superconducting gap. For J \\ c, on the 
other hand, k/T — » as T — » 0. The resulting precipitous decline in the anisotropy ratio k c / K a 
at low temperature rules out a gap structure with line nodes running along the c-axis, such as the 
d-wave state favoured for CeCoIns , and instead points to a hybrid gap of E g symmetry, ft therefore 
appears that two distinct superconducting states are realized in the CeMIns family. 

PACS numbers: 74.70.Tx, 74.20.Rp, 74.25.Fy 



The discovery of magnetically-mediated superconduc- 
tivity in the heavy-fermion material Celn3 [l| has at- 
tracted considerable attention as a possible archetype for 
unconventional pairing in a variety of superconductors. 
However, the fact that the superconducting state in this 
material only exists under pressure makes it difficult to 
know its actual pairing state. Fortunately, the closely re- 
lated family of CeMIn5 compounds offers an ideal test- 
ing ground for investigating the role that dimensionality, 
magnetic order and fluctuations play in determining the 
strength and symmetry of the superconducting state, as 
two members of the family show superconducting order 
at ambient pressure (M = Co, Ir) and the third shows 
antiferromagnetic order [M = Rh) 0, I n CeCoIns, 
the observation of a four-fold anisotropy in the thermal 
conductivity [4j and specific heat 

5] 

on rotation of a mag- 
netic field in the basal tetragonal plane points to a d-wave 
gap (presumably of d x i_ y i symmetry |6j). A number of 
theoretical models propose a d x 2_ y 2 state 0>H,I3' anal- 
ogous to that realized in cuprate superconductors. Since 
the calculated band structure [Hj and measured Fermi 
surface 0] of CeCoIns and Celrlns are very similar, and 
properties like the specific heat [12( and the NQR relax- 
ation rate 0, 0] exhibit the same temperature depen- 
dence, it has generally been assumed that the two su- 
perconductors have the same pairing state, even though 
their transition temperature T c differs by a factor of 6. 
However, because recent evidence suggests that the phase 
diagram of CeMIns may contain more than one super- 
conducting state [lJa, , it has now become crucial to 
pin down the pairing state of Celrlns. 

One of the most conclusive ways to determine the pair- 
ing symmetry of a superconductor is to map out its gap 
structure. A powerful approach to probe the gap struc- 
ture and locate the position of nodes around the Fermi 
surface is to measure quasiparticle heat transport as a 
function of direction, at very low temperature. As an ex- 
ample, the anisotropy of heat transport played a decisive 



role in elucidating the pairing symmetry of the hexag- 
onal heavy-fermion superconductor UPt3 (T c — 0.5 K) 
0,0. In this Letter, we report a study of heat trans- 
port in Celrlns (T c = 0.4 K) down to T c /8 for current 
directions parallel and perpendicular to the tetragonal 
axis of a single crystal. It reveals a dramatic anisotropy 
as T — > 0, whereby low-energy nodal quasiparticles carry 
heat well in the basal plane but poorly, if at all, along 
the c-axis. This is inconsistent with the <i-wave states 
proposed for CeCoIns, characterized by line nodes run- 
ning along the c-axis. In fact, it eliminates all allowed 
(spin singlet) pairing symmetries but one, the (1, i) state 
of the Eg representation. This state has a hybrid gap 
structure, with a line node in the basal plane and point 
nodes in the c-direction, of the kind also found in UPt3 



E3 



Single crystals of Celrlns were grown by the self-flux 
method Q- Two samples were cut into parallelepipeds 
with dimensions ~ 4.5 x 0.14 x 0.045 mm 3 (for J || a) and 
~lx 0.15 x 0.086 mm 3 (for J || c). Their exceptionally 
low residual resistivity (at T — > and H — * 0) attests 
to their very high purity: po a (flo c ) — 0.2 (0.5) fiVL cm. 
The bulk transition temperature is T c = 0.38 ± 0.02 K 
and the upper critical field H C 2 = 0.49 T for H || c. The 
thermal conductivity was measured in a dilution refriger- 
ator using a standard four- wire steady-state method with 
two P1UO2 chip thermometers calibrated in situ against 
a reference Ge thermometer. The same indium contacts 
were used for electrical resistivity and thermal conduc- 
tivity. Their typical resistance at low temperature was 
~ 5 mil. Note that the contribution of phonons to the 
thermal transport is entirely negligible below 1 K. 

Normal state. The thermal conductivity of Celrlns is 
plotted in Fig. ^ as k/T vs T, for a current perpendic- 
ular (J || a) and parallel (J || c) to the c axis. We first 
concentrate on the normal state, where the electrical re- 
sistivity p(T) was found to satisfy the Wiedemann-Franz 
law to better than 1 %, as T — > 0: kn/T = Lq/pq, where 
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FIG. 1: Thermal conductivity n of Celrlns, plotted as k/T 
vs T, for a heat current perpendicular (top) and parallel (bot- 
tom) to the c-axis, in the superconducting state (SC; H = 0) 
and the normal state (N; H — 0.5 T). The normal state data 
are fitted to the Fermi-liquid dependence kjv/T = l/(a + bT 2 ) 
(red line). A linear fit to the superconducting state data for 
J || a (below T c /2) is also shown (blue line). 
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This shows that our measurements do 



j o - —( — ) ■ 
not suffer from electron-phonon decoupling (see discus- 
sion in |l9l Epjp. ktv exhibits the temperature depen- 
dence of a Fermi liquid, k n (T)/T = I /{a + oT 2 ), with a 
= 8.5 (19.6) K 2 cm/W and b = 36 (90) cm/W, for J || a 
(J || c). The fact that both samples have the same (ther- 
mal) resistivity ratio, namely k/T(T -> 0)/k/T(0.6 K) = 
2.4 (2.6) for J || a (J || c), shows that they have the same 
level of impurity scattering. 

Superconducting state. Given this conventional normal 
state behaviour, let us turn to the superconducting state, 
starting with J || a. In the top panel of Fig. ^ n a /T is 
seen to be roughly linear in T, extrapolating to a finite 
residual linear term as T — > 0, as found by Movshovich 
et al. 0| • This is the behaviour expected of a supercon- 
ductor with line nodes, i.e. with a densit y o f states that 
grows linearly with energy {N{E) ex E) [U |H H3| . In 
particular, theory shows that, in the T — > limit, k/T 
reaches a universal value given by [22I |24| : 
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TABLE I: Even-parity (spin-sin glet) pair states in a tetrago- 
nal crystal with point group Da,h |27||- (V = vertical line node, 
H = horizontal line node.) 
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where 7jv is the linear term in the normal state specific 
heat, vf is the Fermi velocity, Ap is the gap maximum, 



fi is the slope of the gap at the node, and a is a con- 
stant of order unity whose value depends on the particu- 
lar gap structure |22|. kq/T is called "universal" because 
it does not depend on impurity concentration and can 
therefore be used to measure the magnitude of the gap. 
Experimentally, universal conduction has been observed 
in high-X'c cuprates |25| and in the spin-triplet supercon- 
ductor Sr2RuC>4 [2^, and Eq. 1 works well in both cases. 

Let us now apply Eq. 1 to Celrln 5 . The allowed order 
parameter representations in tetragonal symmetry |27l | 
are listed in Table I (for singlet pairing). Two line node 
topologies are possible: vertical line nodes (where the 
Fermi surface cuts a vertical plane, e.g. x = 0), such as 
in the two d-wave states (d x 2_ y 2 in B\ g or d xy in B2 g ), 
and a horizontal line node (where the Fermi surface cuts 
the z = plane), such as in the hybrid gap of the E g (1, i) 
state. The simplest gap functions are A = A o cos20 and 
A = 2Aocos0sin0e 4 ^, for d-wave and E g (l,i) states, re- 
spectively. The corresponding nodal structures are illus- 
trated in Fig. [21 Let us apply Eq. 1 to such a hybrid 
gap function, for which a = 3/2 and /1 = /iiinc = 2 [24j . 
Using the known values of jn (7300 J K -2 m -3 H3), vf 
(2x 10 4 m/s, in the basal plane |ll)h and A (2.5 fcsT c 
[13), Eq. 1 yields n 0a /T = 28 mW/K 2 cm. From Fig. □ 
the experimental value is K^ a /T ~ 20 mW/K 2 cm, in 
remarkable agreement with the theoretical estimate. If 
we model the hybrid gap structure in terms of a single 
effective slope of the gap at the line node averaged over 
the various sheets of the Fermi surface of Celrlns, that 
parameter comes out to be /z e ff = 2.8. While this quanti- 
tative agreement with theory is compelling confirmation 
for the presence of a line node in the gap of Celrlns, it 
actually says little about its location. Indeed, the corre- 
sponding estimate for a o!-wave gap gives a similar value 
for KQa/T. The diagnostic power of thermal conductivity 
in determining the detailed topology of the gap comes 
from its directional character, accessed by sending the 
current in distinct high-symmetry directions of the crys- 
tal. This was not done in the previous heat transport 
study and, to the best of our knowledge, no direc- 
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FIG. 2: Typical gap structures in tetragonal symmetry, 
drawn on simple Fermi surfaces. Left : d x 2_ y 2 gap on a 
(warped) cylindrical Fermi surface, with four vertical line 
nodes running along the c-axis. Right : hybrid gap on a spher- 
ical Fermi surface, with one horizontal line node in the basal 
plane and two point nodes along the c-axis. Quasiparticle 
heat conduction at low temperature is entirely governed by 
nodal topology. 

tional measurement of the gap has been reported so far 
for Celrlns. 

Anisotropy. As seen in Fig. 1, applying the current 
along the c axis reveals a qualitatively different limiting 
behaviour, whereby n c /T — > as T — > 0. Simple T 2 or 
T 3 extrapolations yield kqc/T values no greater than 1- 
2 mW/K 2 cm, an order of magnitude smaller than Ko a /T. 
Fig.[3]shows the anisotropy ratio, k c / n a , as a function of 
temperature, in both normal and superconducting states. 
In the normal state, n c / n a is virtually independent of 
temperature, with K a / k c ~ 2.5. The anisotropy in p(T) 
is similarly constant, even well beyond the Fermi- liquid 
T 2 regime, with p c / p a ~ 2.7 between 1.2 and 8 K. This 
simply reflects the anisotropy of the Fermi velocity (or 
mass tensor). A T-independent normal-state anisotropy 
was also found in UPt 3 \mUM 

The superconducting state anisotropy is strikingly 
different, a difference that can only come from gap 
anisotropy. Two distinct features are manifest: 1) a slight 
increase immediately below T c and 2) a precipitous drop 
below T ~ T c /3. These two features combine to produce 
a broad peak centered at T ~ T c /2. We attribute the first 
feature to an anisotropic suppression of inelastic scatter- 
ing, brought about as electrons pair up (anisotropically) 
and cease to participate in the electron-electron scatter- 
ing responsible for the bT 2 term in kn/T. 

The second feature is directly diagnostic of the nodal 
structure, as it comes from low-energy quasiparticles. 
The factor of ~ 3 drop in k c / K a between T c /3 and T c /8 
clearly extrapolates to a very small value as T — > 0. This 
reveals a qualitative a-c anisotropy in the average velocity 
of thermally excited nodal quasiparticles. In other words, 
those fc-states responsible for c-axis conduction in the 
normal state appear to be much more strongly gapped. 
This excludes any nodal structure for which the line nodes 
are along the c-axis, irrespective of the shape of the Fermi 
surface. Indeed, such vertical line nodes would simply 
reproduce the underlying anisotropy of vp, and K c /« a 
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FIG. 3: Temperature dependence of the anisotropy ratio 
« c / n a of Celrlns, in the normal state (N) and in the super- 
conducting state (SC). The dotted line is the ratio of the two 
fit lines (to the normal state data) displayed in Fig. 1 and the 
solid line is a linear fit to the superconducting state data be- 
low T c /3. The precipitous drop at low temperature reflects a 
strongly anisotropic gap whose nodal structure is inconsistent 
with vertical line nodes (running along the c-axis). The small 
peak below T c is due to inelastic scattering (see text). Inset: 
calculated anisotropy (normalized at T c ) for the gaps shown 
in Fig. |21 namely d-wave |2^| and hybrid |2l|. 

would basically mimic the normal state anisotropy. This 
expectation, confirmed by calculations [2^ , is illustrated 
in the inset of Fig. 3 (horizontal blue line). A modu- 
lation of the gap maximum along the c-axis, whereby 
Ao = Ao {&), can produce some additional anisotropy 
in the superconducting state, but this is typically mod- 
est and weakly T dependent j^. More importantly, it 
would never bring to nearly zero for J || c (6 = 0) the 
residual linear term present for J || a (9 = tt/2). 

By excluding vertical line nodes in the gap of Celrlns, 
our study eliminates all allowed representations for the 
order parameter, except one: the two-component E g 
representation (see Table I). In particular, both d-wave 



states are ruled out: cL 



and d xy , respectively in B 



is 



and E>2 g symmetry. Of the three states allowed in the 
Eg representation, only the (1, i) state is generically free 
of vertical line nodes. Its typical [x + iy)z dependence 
produces a hybrid gap, which possesses, in addition to 
the line node in the basal plane (z = 0), point nodes 
along the z \\ c direction, at x = y = (see Fig. 2). 
Note that this state breaks time-reversal symmetry, and 
will therefore spontaneously generate an internal mag- 
netic moment around impurities. /iSR measurements on 
Celrln 5 have not detected such moments 0, possibly 
because the associated fields are too small in these high 
purity samples, as also found in high quality crystals of 
UPt3 30], a superconductor for which the accumulated 
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evidence points overwhelmingly to a ground state with 
broken time-reversal symmetry |18| . 

We now consider whether our data is compatible with 
another special feature of the E g state: the c-axis 

point nodes of its hybrid gap. These are linear point 
nodes, i.e. A(6>) oc 9, such that N(E) oc E 2 , which 
implies that Ko/T in the c direction is not universal. 
Theory shows that kq c /T should be smaller than Ko a /T 
by a factor ~ y/HT/Ao x (i>tine/ ftpomt ( m the unitary 
limit) 24], where T is the impurity scattering rate and 
Mline (Mpoint) is the slope of the gap at the line (point) 
node. The predicted anisotropy for HY = 0.1 ksT c and 
Mine = Mpoint = 2 is shown in the inset of Fig. 3 (red 
line; taken from [!]]])• While the detailed temperature 
dependence of the data on this multi-band material is 
not expected to be captured by the simple model of a 
single spherical Fermi surface, it is nevertheless mean- 
ingful to look at the T — > lim it, gov erned entirely by 
the slope of the gap at the nodes |22j,|24|. The calculation 
shown in Fig. 3 yields a residual anisotropy K c /na that 
is 20 % of the normal state anisotropy at T — * 0. This is 
compatible with the data, where the lowest point only re- 
stricts the residual anisotropy ratio to be less than 40 % 
of its normal-state value. Even a sizable increase in T 
would not necessarily invalidate this compatibility since 
it could easily be compensated by a reduction in the gap 
parameter /i po int- To summarize, while Ko a /T is incon- 
sistent with an s-wave gap {A\ g in Table I) and kq c /T is 
inconsistent with a d-wave gap (or any gap with a verti- 
cal line node), both Koa/T and k,q c /T are quantitatively 
consistent with a hybrid gap. 

We note that the anisot ropy of heat conduction was 
also measured in CeCoIns |2fJ, but the presence of un- 
paired electrons in that material produces an unexpect- 
edly large and isotropic residual linear term which totally 
masks any anisotropy that might come from the coexist- 
ing nodal quasiparticles. It should be emphasized that 
the lack of a sizable residual linear term in the c-axis 
data reported here rules out the possibility of such un- 
condensed electrons in Celrln5. 

In conclusion, the in-plane thermal conductivity n a of 
Celrlns measured down to T c /8 confirms unambiguously 
the presence of line nodes in the superconducting gap. 
The c-axis conductivity n c reveals a profound anisotropy 
as T — > 0, which rules out the possibility that these line 
nodes are vertical (along the c-axis). This eliminates 
all but one of the pairing states allowed in D^h symme- 
try, including the d-wave state proposed for the closely 
related compound CeCoIns. This leaves as sole candi- 
date for Celrlns the state of the E g representation, 
also a prime candidate for the superconductor UPt3 . 
The T — > value of k/T in both high-symmetry direc- 
tions is in good quantitative agreement with theoreti- 
cal calculations for this state. This therefore points to 
superconducting order parameters of different symmetry 
in the two isostructural members of the CeMIns fam- 



ily of nearly magnetic heavy-fermion metals. It will be 
interesting to examine how these differences might arise 
from the respective magnetic fluctuation spectra. Di- 
rect experimental confirmation of the presence of c-axis 
point nodes and broken time-reversal symmetry, both im- 
plied by the E g (l,i) state, is called for. In principle, 
both should be revealed by doping with impurities. We 
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